Multimodal transition and stochastic antiresonance in squid giant axons 
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The experimental data of N. Takahashi, Y. Hanyu, T. Musha, R. Kubo, and G. Matsumoto, 
Physica D 43, 318 (1990), on the response of squid giant axons stimulated by periodic sequence of 
short current pulses is interpreted within the Hodgkin-Huxley model. The minimum of the firing 
rate as a function of the stimulus amplitude Iq in the high-frequency regime is due to the multimodal 
transition. Below this singular point only odd multiples of the driving period remain and the system 
is highly sensitive to noise. The coefficient of variation has a maximum and the firing rate has a 
minimum as a function of the noise intensity which is an indication of the stochastic coherence 
antiresonance. The model calculations reproduce the frequency of occurrence of the most common 
modes in the vicinity of the transition. A linear relation of output frequency vs. Iq for above the 
transition is also confirmed. 

PACS numbers: 87.19.11,87.19.1n,87.19.1c 



The Hodgkin-Huxley (HH) model^ is a prototypical 
resonant neuron with the main resonant frequency typi- 
cally of order 40 to 60 Hz. Its output interspike intervals 
(ISI) can be classified in terms of integer multiples of 
the driving period. The multimodality is revealed when 
the HH neuron is stimulated by noisy inputs, such as 
additive noise|2, 3], random synaptic inputsQ, Q, Q or 
channel noise [6i]. Such ISI histograms are encountered 
frequently in periodically forced sensory neurons. An 
explanation in terms of a two-state system with noise 
was put forward by Longtin et al.0- The multimodal 
character is manifest also in a deterministic HH model 
near excitation threshold [8,9] and in regimes of irregular 
response between mode- locked states [9]. It was shown 
rece ntly that also the parity of ISI plays a significant 
role[10]. Even (odd) modes dominate in the vicinity of 
even (odd) mode-locked states, respectively. The most 
significant manifestation of this effect is the multimodal 
odd-all transition between states 3:1 and 2:1 [10], where 
the coefficient of variation (CV) has a maximum and the 
firing rate has a minimum. The notation p:q means p 
output spikes for every q input current pulses. Below 
this singularity only odd multiples of the input period 
exist and above it harmonics of both parities participate 
in the response. The transition may be crossed by vary- 
ing either the stimulus amplitude or the input period. 
The minimum of the firing rate occurs shghtly above the 
transition. 

In earlier experiments in giant axons of squid stimu- 
lated periodically by a train of short rectangular current 
pulses the firing rate, defined as the ratio of the out- 
put and input frequency fo/fii had a well pronounced 
minimum as a function of the interval between adjacent 
pulses [nj or the stimulus amplitude [l2l|. Even modes 
were absent below the minimum[T2]. This effect occurred 
near the excitation threshold, between states 3:1 and 2:1. 
Another interesting result was the continuous relation be- 
tween the firing rate and the stimulus amplitude. This set 
of experimental and theoretical results deserves a more 



detailed comparison. 

The theory can be tested also by considering a pe- 
riodic drive in the presence of noise. Noisy biological 
systems 0, [3, [1343, including the HH neuron, are known 
to exhibit stochastic resonance (SR). This phenomenon is 
mainly, though not exclusively, characterized by a maxi- 
mum of the signal to noise ratio as a function of the noise 
intensity. Another effect associated with the presence of 
noise is the decrease of the firing threshold and the coher- 
ence resonance [TtI [18| . where the minimum variability of 
the output signal, expressed by CV in absence of a de- 
terministic drive, is achieved at some intermediate noise 
strength. Recently it was found experimentally [l3, [13 
theoretically [21I, [22| that small amplitude noise may de- 
crease the firing rate or even turn it off. The nonlinear 
system in the vicinity of the multimodal transition is a 
natural candidate for finding interesting effects due to 
noise since the trajectories of different modes are very 
close in parameter space. In the following we compare 
experimental data to theoretical results for the determin- 
istic case and calculate the sytem's response to a periodic 
drive with additive Gaussian noise. 

In the experiment of Takahashi et al. [12j the squid axon 
was stimulated by periodic train of rectangular current 
steps of width 0.6ms. Fig. [T] shows the experimentally 
obtained firing rate as a function of stimulus amphtude 
scaled by the minimum current threshold It- On the 
left side of the minimum only odd modes were recorded. 
Even modes were present at the minimum point, with the 
6 : 1 mode occurring more frequently than the 4 : 1 com- 
ponent, and 2 : 1 entirely absent. This is consistent with 
calculation results where even modes disappear be- 
fore reaching the multimodal transition (which is shghtly 
below the minimum of the firing rate), with the low order 
modes vanishing first, beginning with mode 2:1. 

We try to reproduce this type of dependence using 
the HH model with the classic parameter set and rate 
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FIG. 1: The average firing rate, To/Ti, as a function of the 
stimulus amplitude / from the work of Takahashi et al.p!^. 
It is the minimum current threshold obtained in the range 
Ti — 2.5ms to Ti — 6.5ms. The measurements were carried 
out at Ti — 3.8ms. 



constants [ll, 



C — = -iNa - Ik - II + lappi 



(1) 



where In a-, Ik-, Ili lapp, are the sodium, potassium, leak, 
and external current, respectively. C = 1/iF/cm is the 
membrane capacitance. The input current is a periodic 
set of rectangular steps of width 0.6ms and height Iq. 
Equations are integrated within the fourth order Runge- 
Kutta scheme with a time step of 0.001ms. The data 
points are obtained from runs of 400 s, discarding the 
initial 4 s. The dependence of the firing rate on the stim- 
ulus amplitude is shown in Fig. [2l where Ti = 7ms. 




12 14 16 18 20 22 24 
Iq (|lA/cm^) 

FIG. 2: The calculated average firing rate at Ti = 7ms with- 
out noise. 

The similarity to experiment is striking. Although the 
calculated minimum occurs at almost twice the experi- 
mental Ti, the other time scales such as the refractory 
period and the time span of the bifurcation diagram dif- 
fer by a similar factor. The entire dynamics of the axon 
from the study of Takahashi et al. is significantly faster 
than that of Hodgkin and Huxley. This difference of time 
scales is not unusual. Long time ago Best [23,] noted that 



TABLE I: Frequency of occurrence of the six lowest modes at 
the minimum of the firing rate. The upper row is based on Fig. 
13e from the experimental data of Takahashi et al.fl^. The 
bottom row is the result of calculations, assuming Ti — 7ms 
and lo = 18//A/cm^ (see Fig. 
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the axon used by Hodgkin and Huxley was of poor qual- 
ity and in later studies significantly higher conductivities 
were obtained. Paydarfar et al.[19| in their recent study 
recorded firing periods in the range between 7 and 16 
ms. The overall dynamics of Figs. [T] and [2] agrees very 
well, including the location and depth of the local min- 
ima. We verified that the form of Fig. |2]was unchanged 
for pulse widths between and 1ms after dividing the 
current amplitude by f^' I{t)dt. 

Fig. [3] shows the response diagram in the high- 
frequency hmit. The dotted hue separates the monos- 
table firing solution from the silent state and bistable 
areas where the limit cycle coexists with a fixed point so- 
lution. Boundaries of bistability were determined using a 
continuation method starting from a region with a single 
solution. 
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FIG. 3: The bifurcation diagram in the Ti-Io plane, showing 
the main mode-locked states in the model without noise. The 
unmarked intrusion in the upper left corner is the 5:1 state. 
The bottom part of the figure is occupied by the silent state. 
In the firing part of the diagram there are two solutions below 
the dotted line. Here the limit cycle coexists with the fixed 
point. Full squares show the location of the minima of the 
firing rate. The borders of states below Ti — 4.5ms are shown 
in an approximate form. The detailed picture is less regular 
and somewhat more complex. 

The experimental local maximum on the plateau 
fo/fi = 0.4 is due to the state 10100, where modes 
2:1 and 3:1 alternate. The other local maximum at 
f^/f^ = 0.429 with tendency to lock into the (10)^100 
was also reproduced. Fig. [U shows the relative frequency 
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of participation of the most common modes on a logarith- 
mic scale. Higher-order modes appear more frequently 
near the minimum of the firing rate. Experimental and 
calculated ISI histograms are compared in Table HI The 
overall agreement is quite remarkable. Also the calcu- 
lated evolution of individual modes as a function of Iq is 
close to measured values. In experiment the probabihty 
of appearance of mode 4:1 between lo/It = 1-2 and 1.3 
remains in the range 0.06 to 0.08, which agrees well with 
Fig. Hfor Iq between 18/iA/cm^ and 22/iA/cm^ The 
pubhshed experimental runs[12'| contain 80 to 100 output 
spikes for selected data points. On the basis of this data 
set we can conclude that the frequency of participation 
of the low order modes is approximately reproduced in 
simulations. Above the multimodal transition the exper- 
imental firing rate near the threshold rises linearly with 
the stimulus amplitude, see Fig. |5l The dependence of 
fo/ fi vs. of Iq is well reproduced in Fig. [6l except in the 
vicinity of the 2:1 plateau, where an addition of a small 
amount of noise would improve the fit. 
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FIG. 4: The relative frequency of occurrence of low-order even 
and odd modes for the parameter set of Fig. [2] The vertical 
line marks the position of the minimum of the firing rate. 

We now consider the model with a Gaussian white 
noise: 

= -iNa -Ik-Il+ lapp + C^t), (2) 

where < ^i{t) >= 0, < ^(t)^(t') >= 2D5{t - t'), and 
D is expressed in mV^/ms. The HH equations are in- 
tegrated using the second-order stochastic Runge-Kutta 
algorithm [24] . The simulations are carried out with the 
time step of 0.01ms and are run for 400 s, discarding the 
initial 40 s. 

There is a tendency to assume that biological systems, 
including neurons, should always be treated as noisy sys- 




FIG. 5: The linear relation of the firing rate vs. the stimulus 
amplitude above the multimodal transition point at Ti — 4ms. 
These are experimental results of Takahashi et al.ll3]. 
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FIG. 6: Calculated average firing rate vs. stimulus amplitude 
above the multimodal transition for three values of Ti. The 
current pulse width is 0.6ms. 



tems. While the neuron is sensitive to noise it is not ob- 
vious that single neuron dynamics should always include 
stochastic terms. Fig. [T] shows the quick disappearance 
of the fo/fi = 0.4 plateau in Fig. [2]with increasing noise. 
Comparing with the experimental data in Fig. [T]we con- 
clude that calculations reproduce experimental data only 
for D < 10~^. Certainly more experiments are needed to 
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FIG. 7: Sensitivity of the fo/fi = 0.4 plateau from Fig. [2] to 
noise. 
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understand the role of noise in neurons. 

Fig. [8] presents the firing rate as a function of D for 
three parameter sets from the 3 : 1 plateau of Fig. [2l For 
small noise fo/ fi drops quickly below 1/3 over an en- 
tire plateau, with the biggest drop near the edges. This 
behavior should be contrasted with the resonant regime 
where the central part of each plateau maintains phase 
locking over much larger range of noise intensities and 
Z) ~ 1 is needed to lower the firing rate of an entire 
plateau below the Z) = valuefo']. Another difference is 
the direction of frequency changes at the plateau edges. 
In the resonant state the frequency below (above) the 
plateau midpoint is lowered (increased), respectively [9]. 
In the antiresonant limit the entire plateau is unstable to 
even a small noise which slows down the system consid- 
erably. 
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FIG. 8: The firing rate vs. the noise intensity. The middle 
curve was obtained for Iq — 15/xA/cm^. Here Ti — 7ms. At 
D = all three curves start in the 3:1 mode. 

CV as a function of D has a maximum for the same 
parameter set, see Fig. [9l The increased variability is 
associated with increased participation of higher-order 
modes and may be called a stochastic coherence antires- 
onance. A maximum of CV was found earlier in a leaky 
integrate-and-fire model with an absolute refractory pe- 
riod for suprathreshold base current [25j. A small local 
maximum of CV at intermediate noise level was also 
found by Luccioh et al.[5] in a HH model driven by a dc 
current in a bistable regime, where the neuron was stim- 
ulated by a large number of stochastic inhibitory and 
excitatory postsynaptic potentials. It was pointed out 
that the stochastic antiresonance may exist in regions of 
bistability[22j, when the stable limit cycle coexists with 
other attr actors. This typically occurs in the vicinity of 
a bifurcation when the value of the bifurcation parame- 
ter shghtly exceeds the critical value. In the HH model 
near the multimodal transition there are many compet- 
ing limit cycles. Noise enhances trajectory switching and 
may even stop the firing entirely. A decrease of the fir- 
ing rate and an increase of incoherence may occur along 
much of the excitation threshold, where the deterministic 
system is bistable or responds irregularly[9j. 
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FIG. 9: The maximum of the coefficient of variation as a 
function of the noise intensity is a property of the stochastic 
coherence antiresonance. The maximum of CV and the min- 
imum of the firing rate occur at different noise levels. The 
irregularity of the Jo = 14/xA/cm^ curve is a consequence of 
proximity to the excitation threshold. 



In conclusion, numerical solutions of the determinis- 
tic HH equations show that the minimum of the firing 
rate observed by Takahashi et al.[12j is due to the multi- 
modal transition[10j. The statistics of the experimental 
spike trains confirm that below the transition only odd 
modes remain. Even modes are present at the minimum 
of fo/ fi, in agreement with theoretical calculations 0. 
The calculated frequencies of occurrence of the most com- 
mon modes are close to experimental values. Also the 
location of the minimum of fo/ fi in the vicinity of the 
3:1 state is consistent with the simulations. The linear 
rise of the output frequency as a function of the stim- 
ulus strength above the multimodal transition was also 
confirmed. The excitation threshold in the antiresonant 
limit is higher by about a factor of two compared to the 
resonant regime. The rise of threshold for frequencies of 
current pulses exceeding the resonant frequency was ob- 
served experimentally by Kaplan et al. [26j . Further sup- 
port for the significance of the parity of the modes comes 
from the experiment of Paydarfar et al. [l9| , who found 
that the quiescent periods between highly regular bursts 
were always equal to even multiples of the resonant pe- 
riod. An ISI histogram with odd modes was obtained by 
Racicot and Longtin0 in a chaotically forced FitzHugh- 
Nagumo (FHN) model. FHN equations are often used as 
a substitute for the full HH model. It would therefore 
be useful to investigate whether the main features of the 
odd-all ISI transition are reproduced in the FHN model 
with a deterministic and stochastic drive. 

Perturbing the system with noise changes significantly 
the / vs. Iq dependence. The local minima of this curve 
disappear already for D ~ lO""^. In the regime below the 
multimodal transition the 3:1 plateau disappears rapidly 
for very small noise. The firing rate has a minimum and 
CV has a maximum as a function of the noise inten- 
sity. These predictions are expected to be valid for short 
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stimuli of different shapes and can be tested experimen- 
tally. The multimodal transition and the accompanying 
stochastic antiresonance are important both for the un- 
derstanding of excitable systems and for potential neu- 
rological applications. Spike annihilation by determin- 
istic signals 1 28] is studied in the context of deep brain 
stimulation [29j, |30|, a therapeutic technique, in which 
synchronicity of certain parts of the brain is reduced by 
brief current pulses. 
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